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Abstract 

We derive expressions for three-body phase space that are explicitly symmetrical 
in the masses of the three particles, by three separate methods. 
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Nature of the problem 

The phase volume in D-dimensional space-time for the decay process p — > 1 + 2 + 3 can 
be reduced to the integral 
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where s = (pi + P2) 2 , t = (p 2 +P3) 2 , u = (j>3 + Pi) 2 are the three Mandelstam variables 
and, in the physical region, the Kibble cubic 
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stays positive. Clearly p is symmetric in the three masses. However, when one eliminates 
one of the Mandelstam variables (or any linear combination) one is left, for even D, with 
an elliptic integral which is not explicitly symmetrical, although it must be so implicitly. 
The aim of this article is to exhibit three routes for getting a satisfyingly symmetrical 
result in the 4D case, where p is nothing but the area of the Dalitz-Kibble plot 
divided by 128ir 3 p 2 . Detailed derivations are not given, for lack of space, but will be 
published elsewhere. 

Route 1: Elliptic Drive 

The complicated expressions found by Almgren || and Bauberger et al 0, obtained by 
finally integrating [s lj2 = (mi ± m 2 ) 2 , s 3i4 = (p =F m, 3 ) 2 below] 
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involve the complete elliptic integrals, K(k), E(k) and U(a 2 ,k), where 
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The main point is that the result is symmetric under 1 «-> 2 interchange but not obviously 
so under 1 «-> 3 and 2 «-> 3 interchanges. The key to converting these expressions into 
symmetric form is to pass to the Jacobian zeta-functions Z via 
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with sin/3 = a/k and to use the addition formula (142.01) of Byrd and Friedman H: 
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After numerous manoeuvres, including identification of the three angles 
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etc. and noting the remarkable connection 

Z(</?i, A;) + Z(v? 2 , ^) + Z(y2 3 , fc) = ^ 2 sine/?! siny2 2 sin</? 3 , 
one obtains the beautifully symmetric form: 
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Route 2: Phase Cut Way 

In the triangular phase plot, pictured in Figure 1, points O refer to configurations where 
the Mandelstam variables assume their minimum values, points P where they attain 
their maximum values and the Q-lines and i?-lines are tangents at P and O. Needing to 
integrate within the cubic outline OPOPOP, one is guaranteed to obtain a symmetric 
answer by finding the area of the hexagon Aq with Q vertices and cutting out the six 
OQP crescents. Now area 

Aq — — (p — mi — m 2 — m 3) 2 (p + m i + m 2 + rn 3 ) 2 — A(m\ + m\ + m 2 ) , 



while a typical crescent has an area 
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that is definitely given by an elliptic function. These crescents are readily evaluated in 
the limit of large p and of small Q = p — (m 1 + m 2 + m 3 ) by analytical means. One finds 
to leading order 

\ / V \ 31 
log 



^-OtPsQsu — 



1 1 

m 3 p 



m2m 2 p 2 



V 

m 2 + m 3 J 2 

3 , oq 2 , on™2_ , a 3 



12(m 2 + m 3 y 



12m'l + 23m 2 m 2 3 + 20m 2 m 3 + Qm 2 + 0(p) , 



AOtPsQsu 

sin 6»i3 



= 2Q 2 



m 1 m 2 - y l m 1 m2m 3 (mi + m 2 + m 3 ) 6 13 



mim 3 



(mi + m 2 ){m 2 + m 3 ) 



etc. 



3 




Figure 1: Principal extremal points of the Dalitz-Kibble plot 



Subtracting all crescents and noting that 9 12 + #23 + #31 = tt/2 we get the full phase space 
area (in the two respective limits), 
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or 
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Higher order terms can be calculated quite systematically in this manner from the original 
Aqpq integrals for the crescents. 

Route 3: Sunset Boulevard 

This method relies on the observation that ||[7],|| phase space is twice the imaginary part 
of the sunset Feynman integral ||, possessing 3 internal lines, according to Cutkosky's 
cutting rules. Thus p = 21mls(p), where in D-dimensions, 
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The Fourier integral leads to an integral of a Bessel function J of the first kind with a 
product of three K functions of the second kind; balanced progress may be achieved via 
the expansion, 
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which terminates for odd D but is asymptotic for even D. Since for example, 
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we are left with a doable integral, that produces a series of hypergeometric functions, 
whose leading term is 
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with a = rrii + m 2 + m 3> w — mim2m 3 . As imaginary parts of hypergeometric functions 
are || 
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we deduce that 
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(p 2 - a 2 ) = Q(2mi + 2m 2 + 2m 3 + Q) , 

and recognize the threshold behaviour found by other routes. This is a particular case of 
the general iV-body threshold behaviour, Pn(p) ~ q( n - 1 ) d / 2 -( n + 1 )/ 2 _ One last comment: 
it is best to input the massless propagators iA c (x|0) = T(D/2 — 1) / ' (Aix D / 2 r D ~ 2 ) at the 
stage of Fourier integration, rather than attempt to take the zero-mass limit of the above 
massive formulae, because that procedure is extremely tricky. 
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